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Abstract 

We consider a supersymmetric SU(4) c xSU(2)z / xSU(2) j r model with a minimal number of 
Higgs multiplets and a single CP- violating phase 5c p in the neutrino flavor mixing matrix. 
The model incorporates the charged fermion masses and quark mixings, and uses type I seesaw 
to explain the solar and atmospheric neutrino oscillations. We employ thermal leptogenesis 
and the observed baryon asymmetry to find the allowed range for Sep and the reactor mixing 
angle #13. For a normal neutrino mass hierarchy, we find two contours in the sin 2 (26 l i3) - 
Sep parameter space, along which the observed baryon asymmetry is reproduced. For both 
contours, a lower bound sin 2 (2#i3) > 0.03 is obtained. The upper bound turns out to be 
sin 2 (2# 13 ) < 0.25 and 0.35, respectively, both of which exceed the current bound from the 
reactor experiment. The CP- violating phase 5c p can be determined from our results if sin 2 (2# 13 ) 
is precisely measured. For the case of inverted neutrino mass hierarchy the predicted baryon 
asymmetry is not compatible with the observed value. 
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The solar and atmospheric neutrino oscillation have established non-zero neutrino masses 
and mixings between different neutrino flavors. Two mass squared differences and two mixing 
angles relevant for the solar and atmospheric neutrinos are measured with good accuracy. 
The neutrino oscillation parameters to be determined in future experiments include the so- 
called reactor angle 6*13, Dirac CP-phase Sep in the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
mixing matrix, and the ordering of the neutrino mass eigenvalues. The T2K collaboration p] 
has recently reported a non-zero value for # 13 at the 2.5a confidence level. Following the T2K 
results, the MINOS collaboration [2] has reported that their data is consistent with the T2K 
results, while disfavoring 13 = at the 89% confidence level. In the not too distant future, 
the accuracy of measurements for masses, mixing angles and CP-phase in the lepton sector 
should be comparable to that in the quark sector. A precise information of quark and lepton 
mass matrices could provide important clues regarding the origin of fermion masses, flavor 
mixings and CP-violations which, most likely, comes from new physics beyond the standard 
model (SM). 

In order to explain the observed neutrino masses and flavor mixings, we need to extend the 
SM. The type I seesaw mechanism [3] is one of the more promising ways not only to incorporate 
the neutrino masses and flavor mixings but to also explain the tiny of neutrino masses naturally. 
A class of supersymmetric (SUSY) grand unified theories (GUT) has attracted much interest in 
this regard. In addition to providing a resolution of the gauge hierarchy problem, the paradigm 
of SUSY grand unification is also supported by the successful unification of the three SM 
gauge couplings at the GUT scale, Mqut — 2 x 10 16 GeV. Among several possibilities, SO(10) 
unification is one of the more compelling ones, with the quark and lepton multiplets of each 
generation unified in a 16 dimensional representation along with a right-handed neutrino. The 
seesaw mechanism is also automatically implemented, being associated with the breaking of 
SO(10) symmetry to the SM gauge group at Mqut, which is fairly close to the desired seesaw 
scale. 

The so-called minimal SUSY SO(10) model [I] with the minimal set of Higgs multiplets 
(10+126) relevant for fermion mass matrices is a natural extension of non-supersymmetric 
SO(10) models considered a long time ago [5]. Because of the unification of quarks and leptons 
in the 16 representation and the minimal set of Higgs multiplets, the fermion Yukawa matrices 
are highly constrained with the quark and lepton mass matrices related to each other. Note 
that the Higgs 10-plet has been used to implement t-b-r Yukawa unification in SO(10) [5J. 
There have been several efforts within the SO(10) framework to simultaneously reproduce the 
observed quark-lepton mass matrix data as well as the neutrino oscillation data [TJ, [8] . It is quite 
interesting that after the data fitting, essentially no free parameter is left and all fermion Yukawa 
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matrices, in particular the neutrino Dirac Yukawa matrix, are unambiguously determined. The 
neutrino Dirac Yukawa matrix allows us to provide concrete predictions for proton lifetime [9] 
and the rate of lepton flavor violations |10j . 

However, the minimal SO(10) model suffers from a serious problem. The observed neutrino 
oscillation data suggest the right-handed neutrino mass scale to be around 10 13 — 10 14 GeV, 
which is a few orders of magnitude below the GUT scale. With fixed Yukawa couplings of 
right-handed neutrinos in the minimal SO(10) model, this intermediate scale is provided by 
the vacuum expectation value (VEV) of the 126 Higgs multiplet. This indicates the existence 
of many exotic states with intermediate mass scale, which significantly alter the running of 
the MSSM gauge couplings. This has been discussed in [TTJ, where it is shown that the gauge 
couplings are not unified any more, and even the SU(2) gauge coupling blows up below the 
Mgut- To solve this problem, we may extend the minimal model or may consider a different 
direction in constructing GUT models [T2] . 

In this paper we consider a supersymmetric SU(4) c xSU(2) £ xSU(2)^ (4-2-2) model with a 
minimal set of Higgs multiplets and which closely resembles the minimal SO(10) model. The 
Higgs multiplets consist of #1,2,2 : (1,2,2), #15,2,2 : (15,2,2), and #10,1,3 : (10,1,3), which 
correspond to the Higgs multiplets in the minimal SO(10) model, namely #1,2,2 C 10-plet Higgs 
and #15,2,2 + #10,1,3 C 126-plet Higgs. 

The superpotential relevant for the fermion mass matrices is given by 

W = VV , /;/}//:.2.2 + Yi'FiFjHuw + Y^F^H^ (1) 

where F^ : (4, 2, 1) and F : (4, 1, 2) denote the matter multiplets in z-th generation (i = 1, 2, 3). 
Assuming appropriate VEVs for the Higgs multiplets, we can parameterize the fermion mass 
matrices as the follows: 

M u = CiMi, 2)2 + Ci 5 M 15 ,2,2 , 

M d = M 1)2)2 + M 15)2)2 , 
M D = Cl Mi, 2 , 2 - 3 Cl 5M 15 ,2,2 , 

M e = Ml, 2 ,2 - 3M 15 ,2,2 , 

Mr = M 10 ,i, 3 . (2) 

Here M U ,M^, are the mass matrices for up-type and down- type quarks, Mp is the neutrino 
Dirac mass matrix, M e is the charged lepton mass matrix, and Mr is right-handed Majorana 
neutrino mass matrix. They are given in terms of the three fundamental matrices Mi i2j2 , -^15,2,2 
and Mlo,i,3 and the complex coefficients c\ and C15. Note that the relations between fermion 
mass matrices are exactly the same as those in the minimal SO(10) model, except for Mr. 
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In the 4-2-2 model, Mr is independent of the other mass matrices, while it is proportional to 
-^15,2,2 i n the minimal SO(10) model. 

Assuming left-right symmetry 3, the MSSM gauge couplings unify as usual at M GUT , which 
we identify with the breaking scale of 4-2-2 down to the MSSM gauge group. Therefore, the 
procedure for fitting the charged fermion mass matrices is the same as in the minimal SO (10) 
model. On the other hand, Mr being independent of the other mass matrices provides us with 
the freedom to fit the neutrino oscillation data. 

Let us count here the number of free parameters used to fit the charged fermion mass 
matrices. Because of left-right symmetry, M lj2 2 and M 15i2i 2 are 3x3 complex symmetric 
matrices. Without loss of generality, we take a basis where Mi 2 2 is rea l an d diagonal, so that 
the number of free parameters in Mi^ and is 3 + 12 = 15. The two complex parameters 

Ci and C15 introduce an additional 4 degrees of freedom, and therefore in total we have 19 free 
parameters. The degrees of freedom of charged fermion mass matrices are decomposed into 
3 + 6 = 9 for the lepton and quark mass eigenvalues, and another 9 for a unitary matrix for the 
quark mixings which consists of 4 parameters in the CKM matrix and 5 diagonal CP-phases. 
Since the 5 CP-phases are not observable in the SM, we drop these degrees of freedom. Thus, 
we have 14 free parameters to fit 13 observables [H]. In the minimal SO(10) model, this single 
free parameter is adjusted to fit the neutrino oscillation data (see [7] for details). 

Through the type I seesaw mechanism [3], the light neutrino mass matrix is given by 

m v = YlM- R x Y D vl = U MNS D V U T MNS , (3) 

where Yd is the neutrino Dirac Yukawa matrix and v u is the VEV of the up- type Higgs doublet. 
For simplicity, we assume here that all Majorana phases are zero, so that the light neutrino mass 
matrix is transformed by the PMNS matrix Upmns to the diagonal mass eigenvalue matrix D u . 
Using (Eq. [3]), we can express the right-handed neutrino mass matrix as 

m r = v 2 u (y d u; ins d- 1 uI ins yS) . (4) 

Recall that in this 4-2-2 model, Yd is fixed by fitting the Dirac fermion masses and mixings 
in the same manner as the minimal SO(10) model. Therefore, once all neutrino oscillation 
parameters have been measured, we can obtain the complete information of the right-handed 
neutrino mass matrix Mr from Eq. (T4j). Since 6*13 and the Dirac CP-phase 5c p in the PMNS 
matrix have not yet been determined, we obtain Mr as a function of #13 and 8c p from the solar 
and atmospheric neutrino oscillation data. 



1 The left-right symmetry requires us to add -Fj-F) if 10,3,1 to Eq. (TTJ) with a Higgs multiplet #10,3,1 : 
(10,3, 1). This term corresponds to type II seesaw [13] once -ffio,3,i develops a non-zero VEV. Since a more 
complicated Higgs sector seems necessary to induce such a VEV, we do not consider type II seesaw in this 
paper. 
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Models with the seesaw mechanism can also account for generating the observed baryon 
asymmetry in the universe [T5] . 

Y B = ^ = 0.87 x 1(T 10 (5) 

s 

via thermal leptogenesis [TH], where Yg is the ratio of the baryon (minus anti-baryon) density 
(fis) to the entropy density (s). The out-of-equilibrium decays of heavy Majorana neutrinos 
in the presence of non-zero CP-violating phase generates a lepton asymmetry Yl in the uni- 
verse, which is partially converted to the baryon asymmetry through (B+L)-violating sphaleron 
transitions [T7J [TS] . The conversion rate is given by [T^] 

Y = 8N f + AN » Y = --Y (6) 
B 22N f + 13N H L 23 L ' K } 

Here we set Nf = 3 and Nh = 2 for the numbers of fermion families Nf and Higgs doublets 
N H as in the minimal SUSY SM (MSSM). 

The baryon asymmetry produced is evaluated by solving the Boltzmann equations with the 
information of neutrino Dirac Yukawa coupling matrix and Mr. Since Yd is fixed and Mr is a 
function of #13 and Sep in our model, the resultant baryon asymmetry is given as a function of 
#13 and Sep- Therefore, leptogenesis constrains the parameters #13 and Sep so as to reproduce 
the observed baryon asymmetry. 

As mentioned above, the data fitting procedure for the realistic charged fermion mass matri- 
ces is the same as in the minimal SO(10) model, and so in our analysis we employ the numerical 
values in Yd found in [7]. In the basis where the charged lepton mass matrix is diagonal, the 
neutrino Dirac Yukawa coupling matrix at the GUT scale is unambiguously determined and 
explicitly given by 

/ -0.000135 - 0.00273?' 0.00113 + 0.01367 0.0339 + 0.0580?: \ 
Y D = 0.00759 + 0.0119z -0.0270 - 0.00419?: -0.272 - 0.175?: , (7) 
\ -0.0280 + 0.00397?: 0.0635 - 0.01197 0.491 - 0.526?' / 

for tan (3 = 45^. Using Yd, we determine Mr from Eq. (j4j) as a function of #13 and Sep- Since 
the absolute mass spectrum of light neutrinos has not yet been determined, we consider two 
cases for it, the normal hierarchical (NH) case and the inverted hierarchical case (IH). For the 
NH case, the mass eigenvalue matrix D u is given by 

D v = diag ^m , \Jml + Amf 2 , \Jml + Amf 2 + Am| 3 ^ , (8) 



2 Although the output for the neutrino oscillation parameters obtained in [7] is 3er away from the current 
neutrino oscillation data |20j . the experimental data for charged fermion mass matrices are nicely fitted. Since 
Yd is determined only by data-fitting the charged fermion mass matrices, we can safely use this Yd data without 
contradicting any of the experimental results. 



1 



-3 -2 

log 10 [mo (ev; 



Figure 1: Heavy neutrino mass spectrum versus mo for 5c p = (left panel), and versus 5c p 
for mo = 10~ 3 eV (right panel). We have set sin 2 (26 l i3) = 0.11, which is the best fit value in 
T2K results pQ for 5c p = in the normal hierarchical case. 
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Figure 2: Heavy neutrino mass spectrum versus mo for 5cp = (left panel), versus 5cp for 
m = 10~ 3 eV (right panel). We have set sin 2 (26 , 13 ) = 0.14, which is the best fit value in T2K 
results [1] for 5c p = in the inverted hierarchical case. 



while for the IH case 



D v = diae 



rrir 



Am 2 12 + Am^ 3 , J ml + Am|, m 



(9) 



where the lightest mass eigenvalue m is a free parameter. In our analysis, we adopt the 
following values for the neutrino oscillation data [20J: 



Amj 2 = 7.6 x 10 -5 eV 



Am 2 3 



2.4 x 10~ 6 eV 2 



sin 2 (2^ 12 ) = 0.8704, sin 2 (2fl 2 3) = 1.0. 



(10) 



Let us first show the mass spectrum of the heavy Majorana neutrinos (mass eigenvalues of 
Mr) as a function of mo or 5cp- Figured] (left panel) shows the mass spectrum Mi{i = 1,2, 3) 
of the heavy Majorana neutrinos for the NH case as a function of mo with 5cp = and 
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sin 2 (2# 13 ) = 0.11 (the best T2K fit value for 5 CP = 0). Since the VEV of if 10 ,i >3 (which 



the perturbative regime. The right panel shows the mass spectrum as a function of Sep for 
mo = 1CT 3 eV and sin 2 (26*13) = 0.11. The corresponding results for the IH case are shown in 
Figure [2j For the IH case, we find Mi < 10 9 GeV for any value of mo, sin (2 $13) and 5cp- 

As shown in Figures [1] and [21 the heavy neutrino masses are hierarchical for both the NH 
and IH cases. The lepton asymmetry in the universe in this case is dominantly produced by the 
lightest heavy neutrino decay, since the asymmetry produced by heavier neutrino decays are 
almost completely washed-out [21] . Thus, we consider the lepton asymmetry produced by only 
the lightest heavy neutrino decay. In addition, there is a lower bound on the lightest heavy 
neutrino mass, Mi > 10 9-10 GeV, to produce the desired amount of baryon asymmetry [22]. 
For the IH case, the lightest heavy neutrino mass is always found to be below this bound, and in 
our numerical analysis we find that the resultant baryon asymmetry is too small in comparison 
to the observed baryon asymmetry. Therefore, in the following, we present our results only for 
the NH case. 

Since our model is supersymmetric, we need to consider the lepton asymmetry generated 
by the decays of both the lightest heavy neutrino and sneutrino. From Figure [IJ the lightest 
heavy neutrino mass is far below Mqut, and so the effective theory for leptogenesis contains the 
MSSM and three light neutrinos, as well as the lightest heavy neutrino superfield. Although the 
complete Boltzmann equations for this system is quite involved (see [23] for complete formulas), 
because of supersymmetry the lepton asymmetry stored in the SM particles is exactly the same 
as that stored in the sparticles [23]. Since the heavy neutrino mass scale is much higher 
than the typical sparticle mass scale ~ TeV, our system is supersymmetric to a very good 
approximation. Among the many decay and scattering processes involved in the Boltzmann 
equations, it is known that the (inverse) decay process of the lightest heavy (s)neutrino plays 
the most important role in determining the resultant baryon asymmetry, while the others 
are negligible in most of the parameter space [21]. Including only the decay process greatly 
simplifies the Boltzmann equations, so that for the heavy neutrino they are exactly the same 
as in the non-supersymmetric case: 



where is the yield (the ratio of the number density to the entropy density s) of the lightest 
heavy neutrino, Y^ q is the yield in thermal equilibrium, the temperature of the universe is 



breaks 4-2-2 down to MSSM) is Mqut, we require m > 10 4 eV in order to keep within 




(11) 
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normalized by the mass of the heavy neutrino z = Mi/T, H(Mi) is the Hubble parameter 
at T = Mi, and Yj Jj is lepton asymmetry stored in the SM particles. The CP-asymmetry 
parameter, ei, is given by [21] 

1 -— £lm [(nyj)?,.] /(M 2 /M 2 ), (12) 



2ir(Y v Y l 



where 



fix) = vW^)+2^- 
\ x J x — 1 



(13) 



and Y v is the neutrino Dirac Yukawa coupling matrix in the basis where both the charged lepton 
matrix and M R are diagonalized. Using Eqs. (TjJ and (J7J), we can obtain Y v as a function of # 13 
and Sep- 

The space-time density of the heavy neutrino decay in thermal equilibrium, 7^ is given by 
where K\ and K2 are the modified Bessel functions, and 

r Nl = { ^^M l (15) 

is the decay width of the heavy neutrino. Then, we solve the Boltzmann equations with the 
boundary conditions Fjvi(O) = Y^(0) and 3rx / (0) = 0. The lepton asymmetry generated by the 
right-handed neutrino decays is converted to the baryon asymmetry via the sphaleron process 
with the rate of Eq. fl6]) and hence, we evaluate the resultant baryon number as 

Yb = ~^ Y L f (oo) x 2, (16) 

where the factor 2 takes into account the baryon number stored in sparticles. 

For various values of sin 2 (2# 13 ) and Sep, we numerically solve the Boltzmann equations. In 
our analysis, we fix m = 10~ 3 eV. For m < 10 -2 eV, Mi is almost independent of tuq, and we 
find that the results for the generated baryon asymmetry are almost the same. Figure shows 
the resultant baryon asymmetries as a function of Sep for three different values of sin 2 (2# 13 ) = 
0.025 (dashed), 0.1 (solid) and 0.3 (dotted), along with the observed value (red horizontal line). 
The allowed region for sin 2 (2^13) and Sep compatible with the observed baryon asymmetry is 
depicted in Figure HI The desired baryon asymmetry is reproduced along the two contours. 
Interestingly, we find the lower bound sin 2 (2#i3) > 0.03. The upper bound of 0.15 on sin 2 (2#i3) 
from the reactor experiment [20] is more severe than what we find via baryon asymmetry. 
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Figure 3: Baryon asymmetry as a function of 5cp for different values of sin 2 (2#i3) = 0.025 
(dashed), 0.1 (solid) and 0.3 (dotted). The red horizontal line shows the observed baryon 
asymmetry. 
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Figure 4: Contours in (sin 2 (2^13) , 5cp)-plane along which the observed baryon asymmetry is 
reproduced. 
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In the above analysis, as previously mentioned we have dropped all Majorana phases and 
kept only the Dirac phase. The general neutrino mixing matrix is given by U — Upmns x 
diag(e 4flil , e"^ 2 , 1) with two Majorana phases. We have also analyzed this more general scenario 
and found that the observed baryon asymmetry is realized in this case even for #13 = 0. 
Therefore, our ansatz that <fti = <p2 = is crucial in obtaining the bound sin 2 (26*13) ^ 0.03. 

In summary, we have considered a supersymmetric SU(4) c xSU(2)lxSU(2)# model with a 
minimal number of Higgs multiplets and a single CP violating parameter 5c p in the neutrino 
flavor mixing matrix. The model has the same structure in the Yukawa couplings for the 
charged fermions as the supersymmetric minimal SO (10) model, so that the neutrino Dirac 
Yukawa coupling matrix is unambiguously determined by fitting the experimental data for 
charged fermion mass matrices. Using the seesaw formula with the neutrino Dirac Yukawa 
coupling matrix, the right-handed Majorana neutrino mass matrix is given as a function of 
9i3 and 5c p- We have employed leptogenesis and the observed baryon asymmetry to identify 
the relation between #13 and 5c p- Only the normal hierarchical case for the light neutrino 
mass spectrum can reproduce the observed baryon asymmetry, and we find a lower bound of 
0.03 on sin 2 (26 ) 13 ). The recent results from T2K and MINOS indicate a non-zero whose 
value should be more precisely known within the next few years. If #13 is sizable, 5c p and the 
pattern of mass hierarchies can be determined in future experiments. Our results therefore 
have interesting implications for future neutrino oscillation experiments. 
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